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ABSTRACT 


This  memorandum  is  concerned  with  the  problem  of  selecting  a  capacity 
value  for  the  buffer  between  two  asynchronous  production  operations.  The 
methodology  study  is  motivated  by  shortcomings  in  present  methods. 

Objectives  for  this  study  are  that  methods  should  be  at  least  as  efficient 
as  stochastic  simulation  and  that  they  should  provide  a  means  of  examining 
both  the  utilization  of  the  buffer  and  the  productivity  of  the  manufacturing 
process. 

Two  methods  were  developed  and  implemented  in  computer  programs.  Both 
methods  use  the  theory  of  Markov  processes.  The  first  method  (GS.BUF) 
calculates  steady-state  probabilities  for  all  states  of  a  simple  production 
system.  For  certain  types  of  processes  the  results  of  this  model  are  exact. 
The  results  offer  a  good  approximation  for  many  processes.  This  model 
calculates  the  system  productivity  explicitly,  providing  an  opportunity 
for  economic  tradeoffs  between  buffer  capacity  and  other  parameters.  The 
second  model  (BUF.CAP)  focuses  on  the  dynamics  of  the  filling  and  emptying 
of  the  buffer,  under  the  assumption  of  statistical  independence  between 
states  of  the  two  operations.  This  model  admits  the  possibility  of  several 
machines  working  in  parallel  at  each  operation. 
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3.  Background 

Over  the  past  year  I  have  been  involved  in  a  methodology  study  concerned 
with  estimation  of  the  capacity  of  production  lines.  This  study  [Ref  a]  has 
produced  a  general  computer  simulation  (TANDEMT)  capable  of  performing 
stochastic  simulation  on  manufacturing  systems  having  quite  general  structures. 
By  contrast,  the  problem  addressed  by  this  memorandum  is  quite  restricted  in 
scope.  One  of  the  objectives  of  the  manufacturing  productivity  study  is  to 
develop  the  means  for  understanding  the  significance  of  equipment  and/or 
procedural  changes  on  the  productivity  of  a  specific  manufacturing  line. 

The  methods  discussed  here  concern  only  the  effect  of  buffer  capacity  in  a 
simple,  asynchronous  system  with  two  operations  in  tandem  and  an  intervening 
buffer. 

4.  Motivation 


Ref  b  mentions  shortcomings  of  procedures  for  sizing  the  buffer  between 
manufacturing  operations.  An  ad  hoc  method  is  proposed  there  for  generating 
a  reasonably  sized  buffer,  suitable  for  use  as  an  input  datum  to  a  stochastic 
simulation  of  a  developmental  production  process.  The  method  is  not  claimed 
to  be  logically  rigorous.  The  proposal  yields  a  single  number,  but  lacks  a 
measure  of  sensitivity  of  the  process  output  (or  buffer  performance)  to 
buffer  size.  A  sound  approach  to  this  problem  is  needed  which  permits  the 
calculation  of  the  advantage  of  increasing  buffer  capacity.  A  proper  method 
should  permit  efficient  tradeoffs  to  be  made  between  buffer  capacity  and 
other  process  parameters.  The  methods  presented  in  this  memorandum  possess 
these  attributes. 

5.  Definitions 


The  simplest  tandem  system  consists  of  two  operations  running 
asynchronously  with  an  intervening  buffer.  In  this  MFR  the  term  "simple 
production  system"  refers  specifically  to  this  kind  of  system.  More 
complex  systems  may  be  viewed  as  arrangements  of  these  simple  systems. 
Additionally,  the  term  productivity  is  used  here  somewhat  restrictively. 
Productivity  is  a  measure  of  the  efficiency  of  a  production  system  to  produce, 
within  imposed  machine  limits.  As  used  here,  productivity  is  defined  as  the 
ratio  of  the  average  quantity  produced,  in  steady  state,  to  the  production  of 
a  perfect  system  having  the  same  machine  rates.  The  machines  comprising  the 
production  system  are  viewed  in  a  general  way.  They  can  consist  of  automatic 
hardware  or  of  humans  with  simple  tools  or  anything  in  between.  A  property 
of  a  machine  is  that  it  fails  or  requires  adjustment  by  a  repairman  at  random 
operating  intervals.  The  mean  time  between  "failures"  of  a  machine  type  is  a 
property  of  this  machine,  abbreviated  MTBF.  Similarly,  the  mean  time  to  repair 
is  denoted  MTTR.  To  simplify  the  analysis  without  undue  loss  of  generality,  it 
is  assumed  that  sufficient  repairmen  are  available  so  that  essentially  no  time 


*  Alternatively,  at  least  m  repairmen  are  present,  where, 
P{m  or  less  repairmen  busy}>0.9. 
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is  spent  by  machines  waiting  in  a  repair  queue.  This  is  what  is  meant  by 
a  well  maintained  system. 


6.  Objectives 

The  principal  objective  of  this  study  is  to  develop  a  method  for  sizing 
the  buffer  in  a  simple  production  system.  The  method  should  be  at  least  as 
numerically  efficient  as  stochastic  simulation.  Further,  the  method  should 
provide  a  means  of  examining  both  the  utilization  of  the  buffer  and  the 
productivity  of  the  manufacturing  process. 

7.  Methodology  Overview 


Before  getting  into  details  of  the  analysis,  we  consider  the  general 
approach  to  the  methods  of  this  memorandum.  Both  of  the  methods  are  based 
on  the  theory  of  Markov  processes.  The  computer  programs  which  implement 
these  methods  are  found  in  the  Annexes.  In  the  first  method,  with  program 
name  GS.BUF,  the  three  components  of  a  simple  production  system  are  viewed 
as  operating  together  to  generate  various  system  states.  The  system  states 
are  defined  in  terms  of  the  admissible  states  of  the  1st  machine  operation, 
of  the  buffer,  and  of  the  2nd  machine  operation.  As  shown  in  Figure  1,  the 
1st  operation  is  considered  a  single  machine  whose  states  are  (a)  under 
repair,  (b)  waiting  for  a  space  in  the  buffer  to  place  a  completed  part, 
and  (c)  operating  on  a  part.  These  substates  are  numerically  coded  as  0, 

1,  and  2  respectively.  The  state  of  the  buffer  is  just  the  number  of  parts 
occupying  it.  The  admissible  buffer  substates  are  integers  from  0  to  the 
buffer  capacity.  The  2nd  operation  is  also  considered  a  single  machine 
whose  states  are  coded  0,  1,  and  2  for  (a)  under  repair,  (b)  waiting  for 
a  part  to  remove  from  the  buffer,  and  (c)  operating  on  a  part. 


1  waiting  on 
buf  condition 


buf .cap 


1  waiting  on 
buf  condition 


2  up 

Figure  1.  States  of  a  Simple  Production  System 


2  up 
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8.  A  state  of  the  system  is  denoted  by  three  integers  each  of  which 
characterizes  one  of  the  consecutive  components  of  this  system.  Thus, 
state  0,  0,  1  indicates  that  operation  1  is  under  repair,  that  the  buffer 
is  empty,  and  that  operation  2  is  waiting  for  a  part  to  operate  upon.  For 
convenience  the  states  are  labeled  with  a  single  index  i.  The  probability 
that  the  system  occupies  the  i  th  state  at  time  t  is  denoted  p-j(t).  The 
linear  equations  which  functionally  relate  d(pj(t))/dt  to  the  various 
probabilities  of  state  occupancy,  for  all  states,  are  called  the  Kolmogorov 
equations.  These  differential  -  difference  equations  can  be  simply  written 
under  the  assumptions  of  exponentially  distributed  --  time  to  fail,  time 

to  repair,  and  machine  operation  (or  service)  times.  (As  shown  later  these 
assumptions  are  not  too  restrictive.) 

9.  In  stochastic  steady  state  all  the  derivatives  are  set  to  zero.  The 
resultant  set  of  linear  algebraic  equations  is  solved  for  the  state 
probabilities  in  steady  state.  Certain  states  are  identified  which 
collectively  represent  interesting  conditions.  For  example,  those  states 
in  which  operation  1  is  waiting  or  operation  2  is  waiting  for  the  buffer, 
etc.  A  condition  of  particular  interest  is:  operation  2  is  either  waiting 
or  down  for  repair.  The  probability  that  this  condition  obtains  is  the 
fraction  of  time  that  the  simple  system  is  nonproductive.  The  l's 
complement  of  this  probability  is,  thus,  system  productivity.  The 
probabilities  of  these  state  conditions  can  be  used  in  tradeoff  analyses 

to  size  the  buffer.  Economic  factors  can  be  invoked  to  determine  the 
value  of  increased  productivity  versus  the  cost  of  additional  buffer 
capacity.  This  is  the  basis  of  the  first  method  for  sizing  the  buffer. 

10.  A  Second  Approach 

Another  approach  to  buffer  sizing  is  considered.  Altho  lacking  in  the 
comprehensiveness  of  the  first  approach,  it  does  provide  an  approximation 
of  the  probability  that  a  specific  number  of  buffer  spaces  would  be  required 
rf  the  first  and  second  machine  operations  were  not  constrained  to  wait  for 
a  buffer  condition.  Unlike  the  first  approach,  this  method  explicitly 
accounts  for  the  possible  existence  of  several  identical  machines  working 
in  parallel  at  each  machine  operation.  The  basis  of  the  2nd  approach  is  to 
consider  the  states  of  the  first  operation  to  be  statistically  independent 
of  the  states  of  the  second.  Independence  is  a  reasonable  assumption  if 
the  buffer  capacity  is  large  and  if  there  are  no  common  causes  of  machine 
failure.  Then,  the  Kolmogorov  equations  for  the  first  operation  and  for 
the  second  are  independent  and  have  the  same  simple  form.  Generally,  an 
operation  has  N  machines  operating  in  parallel  and  possesses  N+l  states. 

In  this  case  the  state  value  of  an  operation  is  just  the  number  of 
machines  operating.  In  a  two-operation  system  with  Nj  machines  in 
operation  1  and  N2  machines  in  operation  2,  there  are  (Nj+1)(N2+1) 
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system  states.  The  method,  implemented  in  BUF.CAP,  starts  with  the  system 
conditionally  in  each  of  its  states  in  turn,  and  directly  solves  the 
differential  -  difference  (Kolmogorov)  equations  to  obtain  the  time- 
dependent  average  (expected)  production  from  each  of  the  operations, 
conditioned  by  the  initial  state.  Because  of  the  independence  between 
operations,  the  actual  numerical  procedure  considers  each  operation  by 
itself  and  obtains  the  N+l  conditional  trajectories  —  time-dependent 
state  probabilities  and  associated  expected  production  --  storing  the 
results.  The  difference  in  expected  production  from  operations  1  and  2 
at  time  t  represents  the  expected  value  of  the  parts  which  would  be 
added  to  the  buffer,  if  product  1  >  product  2,  or  would  be  removed  from 
the  buffer,  if  product  1  <  product  2.  Recall  that  this  expected  difference 
is  conditional  upon  the  initial  state  condition  (IC).  But  because  the  IC's 
are  random  the  expected  production  differences  are  actually  random  variables. 

11.  These  differences  are  calculated  for  all  system  states  at  a  time  large 
enough  to  allow  the  state  probabilities  to  approach  their  steady-state 
values  (about  4*MTTR).  Then  the  expected  differences  are  rank  ordered 
from  smallest  to  largest.  The  probability  that  the  system  initial  state 
would  be  occupied  in  steady  state  is  calculated  (via  a  product  of 
binomial  probabilities).  These  probabilities  are  associated  with  each 

of  the  ordered  production  -  difference  values.  These  probability  densities 
are  accumulated  to  yield  the  cumulative  distribution  function  (c.d.f.)  for 
the  production  difference  expected  under  these  conditions.  The  program 
BUF.CAP  displays  this  c.d.f.  In  a  balanced  production  system  the  expected 
value  of  the  production  difference  from  this  distribution  is  zero.  The 
mean  and  variance  of  the  random  variable  from  this  distribution  are  calcu¬ 
lated  in  BUF.CAP.  The  variance  depends  upon  the  number  of  machines  at  each 
operation',  the  machine  rates,  and  the  MTBF's  and  MTTR's.  To  reduce  the  risk 
of  * ausing  machines  to  wait  for  the  buffer,  the  required  buffer  capacity  is 
set  to  the  range  in  expected  production  difference  plus  4  standard  deviations. 
This  somewhat  arbitrary  assignment  provides  a  reasonably  small  risk  that  the 
calculated  buffer  capacity  would  be  inadequate.  A  measure  of  risk  is 
provided  explicitly  in  this  approach,  but  no  measure  of  productivity  is 
given  here.  However,  using  stochastic  simulation  I  have  noted  that  the 
marginal  change  in  productivity  at  the  calculated  buffer  capacity  is  about 
0.04  percent  per  percent  change  in  capacity. 

12.  General  Conclusions 


The  most  general  conclusions  concerning  the  above  methods  are  presented 
here.  Details  are  presented  in  later  parts  of  this  memorandum.  In  production 
systems  having  operation  times  which  are  exponentially  distributed  and  in 
which  repair  times  are  exponentially  distributed,  a  steady-state  Markov 
model  of  a  simple  production  system  provides  a  satisfactory  approach  to 
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sizing  the  buffer.  A  model  of  this  type  is  implemented  in  the  computer 
program  GS.BUF.  Even  when  the  machine  operating  rates  are  constants  (or 
nearly  so),  GS.BUF  provides  a  good  approach,  providing  that  several 
machines  are  working  in  parallel  at  each  operation.  When  constant-rate 
single-machine  operations  obtain,  stochastic  simulation  is  the  preferred 
approach  to  buffer  sizing.  The  use  of  GS.BUF  permits  one  to  invoke  economic 
factors  in  sizing  the  buffer.  Machine  operating  rate  (or  operating  time 
interval),  repair  time,  and  buffer  capacity  are  all  important  variables 
which  affect  productivity.  If  one  wishes  to  examine  the  (filling  and 
emptying)  dynamics  of  the  buffer  in  detail,  the  method  implemented  in 
BUF.CAP  is  suggested.  This  method  is  also  useful  in  obtaining  a  point 
estimate  of  buffer  capacity  when  a  more  elaborate  economic  tradeoff  is 
not  appropriate. 

13.  Derivation  of  Equations  for  GS.BUF 

The  general  outline  of  the  theory  for  GS.BUF,  provided  in  paragraph  7  ff . , 
will  be  detailed  here.  To  make  the  exposition  simple,  consider  a  simple 
production  system  with  a  buffer  having  only  3  spaces.  The  theory  of  Markov 
processes  can  be  applied  directly  to  the  states  defined  in  Figure  1  providing 
the  operating  times,  time  intervals  between  failure,  and  repair  times  are  all 
exponential  random  variables.  Other  substates  could  be  added  to  accommodate 
other  distributions  of  these  random  variables.  For  the  present,  consider 
only  the  states  given  in  Table  1. 

14.  Transitions  occur  between  these  states.  The  state  numbers  in  Table  1 
appearing  at  the  left  in  sequence  are  used  as  indices  to  designate  the 
probability  of  state  occupancy.  Thus,  pj(t)  refers  to  the  probability  that 
the  system  is  in  state  1  at  time  t.  The  states  to  which  a  particular  state 
may  transition  are  listed  in  the  column  labeled  "transition-to  states." 
Similarly,  the  states  from  which  transfers  may  occur  are  listed  in  the 
column  labeled  "transition- from  states."  Ordinarily,  Markov  processes  may 
be  represented  diagrammatical ly  by  a  graphical  network  with  nodes  as  states 
and  arcs  as  transitions  —  the  state-transition  diagram.  The  rate  parameters 
for  the  transitions  are  affixed  to  the  corresponding  arcs.  Because  of  the 
visual  complexity  of  the  graph  for  this  process,  only  a  partial  state- 
transition  diagram  is  shown  in  Figure  2.  In  this  case  the  first  seven 
states  are  isolated,  and  all  states  connecting  each  of  these  states  are 
shown  separately.  The  rate  parameters  shown  in  Figure  2  --  r,  X,  y  — 

are  indexed  with  a  1  or  2  to  indicate  the  operation  to  which  it  belongs. 

For  example,  for  state  1,  transitions  to  state  2  occur  at  the  "birth" 
rate  Aj,  which  characterizes  operation  1. 
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TABLE  1 


DEFINITION  OF  STATES  OF  A  SIMPLE  PRODUCTION  SYSTEM 


Example  with  buffer  capacity  of  3  spaces. 


State 

4 

State 

Definition* 

Trans it ion- to 

Transition-from 

Number 

Opn  1 

Buf 

Opn  2 

States 

States 

1 

0 

0 

1 

2 

2,  4 

2 

2 

0 

1 

1,  6 

1,  6 

3 

0 

0 

0 

4,  5 

4,  5 

4 

0 

0 

2 

1,  3,  6 

3,  6,  8 

5 

2 

0 

0 

3,  6,  9 

3,  6 

6 

2 

0 

2 

2,  4,  5,  10 

4,  5,  10 

7 

0 

1 

0 

8,  9 

8,  9 

8 

0 

1 

2 

4,  7,  10 

7,  10,  12 

9 

2 

1 

0 

7,  10,  13 

5,  7,  10 

10 

2 

1 

2 

6,  8,  9,  14 

6,  8,  9,  14 

11 

0 

2 

0 

12,  13 

12,  13 

12 

0 

2 

2 

8,  11,  14 

11,  14,  16 

13 

2 

0 

0 

11,  14,  17 

9,  11,  14 

14 

2 

2 

2 

10,  12,  13,  18 

10,  12,  13,  18 

15 

0 

3 

0 

16,  17 

16,  17 

16 

0 

3 

2 

12,  15,  18 

15,  18 

17 

2 

3 

0 

15,  18,  19 

13,  15,  18 

18 

2 

3 

2 

14,  16,  17,  20 

14,  16,  17,  20 

19 

1 

3 

0 

20 

17,  20 

20 

1 

3 

2 

18,  19 

18,  19 

*  For  operations  1  and  2,  the 
following  meanings: 

integers 

in  the  state  defin- 

ition  have  the 

0 

means  "down" 

or  under 

repair,  with  a  part  being  held. 

1  means  waiting  for  a  buffer  condition— to  place  a  part,  for 
operation  1,  and  to  remove  a  part,  for  operation  2. 


2  means  "up"  or  working  on  a  part. 

The  integer  characterizing  the  buffer  state  is  the  number  of  parts  in  the 
buffer. 
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Figure  2.  Partial  State  Transition  Diagram  for  a  Simple  Production  Process 


DRSMC-SAS  (R) 

SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 


Transitions  to  state  1  occur  from  states  2  and  4.  The  transition  from 
state  2  occurs  at  the  "death"  rate  y^,  associated  with  operation  1.  The 
transition  from  state  4  occurs  at  the  operating  rate  r2  of  operation  2. 
The  birth  rate  for  an  operation  is,  of  course,  the  rate  at  which  that 
operation  is  restored  to  operational  condition  given  that  it  is  "down". 
Thus, 


Ai  =  1/MTTR.  ,  i  =  1,2  ,  (1) 

for  operations  1  and  2.  Similarly,  the  death  rate  for  an  operation  is 
the  rate  at  which  a  failure  occurs,  given  that  the  operation  is  "up". 
Thus, 


y.  =  1/MTBF.  ,  i  =  1,2  .  (2) 

The  machine  rates  (  reciprocals  of  mean  service  time)  are  denoted  by  r, 
and  r£»  for  operations  1  and  2  respectively. 

15.  With  the  aid  of  the  state  transition  diagram,  writing  the  Kolmogorov 
equations  is  a  quite  mechanical  task.  For  example,  as  is  customary  in 
deriving  transition  equations,  consider  a  small  time  increment  h.  Then, 
the  probability  that  state  1  is  occupied  at  time  t  +  h  is  given  as 

Plstate  1  at  t+h}  =  P{no  transition  in  h  occurs  from  state  1,  given 
occupation  at  t) 

*P{state  1  at  t)  +  Pltransition  from  state  2,  given  occupation  of 
state  2  at  t} 

*P{state  2  at  t>  +  Pltransition  from  state  4,  given  occupation  of 
state  4  at  t) 

*P{state  4  at  t}. 

Using  the  abbreviated  notation,  this  expression  becomes 

Pj  (t+h)  =  (l-AjhJp^t)  +  y1hp2(t)  +  r2hp4(t).  (3) 

Then, 

[  p1(t+h)-p1(t )  ]/h  =  -Ajp^t)  +  y1P2(t)  +  r2p4(t).  (4) 

Taking  the  limit  as  h  approaches  zero  and  omitting  the  functional  dependence 
upon  t, 

Pi  =  -A.  p1  +  yj  p2  +  r2  p4  .  (5) 
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SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 

So  much  for  conventional  derivations!  This  expression  can  be  obtained  directly 
from  the  state  transition  diagram  by  writing  as  terms  the  probabilities  of  all 
states  which  transition  to  a  particular  state--the  1st  here--on  the  right  with 
their  transition  rates  as  positive  coefficients.  The  rates  of  all  transitions 
from  the  particular  state  are  collected  and  the  negative  of  this  sum  is  the 
coefficient  of  the  particular  state.  The  Kolmogorov  state  equations  for  the 
simple  system  with  3-space  buffer  are  written  compactly  in  Table  2.  The  first 
column  in  this  table  is  the  index  (i)  of  the  left  hand  side  p^ . 

16.  For  a  simple  system  with  a  buffer  capacity  of  3-spaces,  there  are  20 
states.  The  display  of  the  corresponding  20  equations  is  awkward,  but  still 
manageable.  With  increasing  buffer  size,  writing  the  state  equations 
explicitly  is  infeasible.  Fortunately,  this  display  is  not  necessary.  For 
notational  simplification,  let  £  be  a  column  vector  of  the  state  probabilities 
with  ns  (number  of  states)  elements.  Let  B  be  a  square  matrix  of  coefficients 
with  ns  rows.  Then,  the  Kolmogorov  equations  for  a  simple  production  system 
can  be  written,  generally,  as 

£  =  B£  .  (6) 

Because  the  sum  of  the  elements  of  £  is  unity,  B  is  not  of  full  rank.  Thus, 
equation  (6)  is  not  solved  directly.  In  obtaining  a  solution,  however,  it  is 
useful  to  generate  the  elements  (b-jj)  of  B.  Because  B  is  a  stochastic  matrix, 
the  sum  of  each  of  its  column  vectors  is  zero  (reflecting  the  fact  that  each 
arc  in  the  state  transition  diagram  both  leaves  and  enters  a  node).  This  fact 
is  exploited  to  check  the  validity  of  the  equations  actually  solved  in  GS.BUF. 

17.  In  Table  1,  note  that  the  first  2  states  are  waiting  (i.e.,  1)  states  for 
operation  2,  whereas  the  last  2  states  are  waiting  states  for  operation  1.  In 
each  of  the  states  between  the  2nd  and  2nd  to  last,  a  regular  pattern  is 
observed.  For  a  given  buffer  state  two  0  states  are  assigned  operation  1, 
with  operation  2  taking  the  values  0  and  2.  Next,  two  2  states  are  assigned 
for  operation  1,  with  operation  2  again  taking  the  values  0  and  2.  This 
pattern  of  four  system  states  is  followed  for  each  value  of  the  buffer  state. 
Thus,  with  a  buffer  capacity  m,  the  number  of  states  is 

•  ns  =  4(m+l)+4 

or 

ns  =  4(m+2)  .  (7) 

Because  of  the  regularities  in  the  state  definition  noted  above,  there  are 
regularities  in  the  equations,  which  permit  the  elements  of  the  B  matrix  to 
be  written  recursively.  Starting  with  state  11  (11th  row  of  B) , 

b-jj  =  bi_4}  j_4  ,  11  £  i  £  ns-5  , 

1-4  £  j  £  ns  .  (8) 

This  fact  greatly  simplifies  the  process  of  generating  the  elements  of 
the  B  matrix. 
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TABLE  2 


STATE  EQUATIONS  FOR  A  SIMPLE  SYSTEM 
WITH  3- SPACE  BUFFER 


(See  Table  1  for  definition  of  states). 


Derivative  of 

P {State  No.} 

Terms  in  Right-Hand-Side  Sum 

1 

"Alpl  5  ylp2  5  r2P4 

2 

A1P1  *  "^yl+rPp2  5  r2P6 

3 

™*  ( A i+Ag )  P3  3  3  y  1 P5 

4 

^2^3  9  ”(^l+^2+r2^4  9  ^1^6  9  r2^8 

5 

A1p3  ’  “ (A2+y i+rl^p5  5  y2p6 

6 

r^p2  ,  ,  A^Pj.  ,  -(y^+y^rj+r^Pg,  ^2P10 

7 

-(Aj+A2)p7  9  ^2P8  5  ylp9 

8 

A2P7  ’  "^Al+y2+r2^P8  5  ylp10  »  r2P12 

9 

rlP5  »  A1P7  ’  -(A2+yi+ri^ pg  »  y2p10 

10 

rlP6  '  A1P8  ’  A2P9  5  “(yi+y2+ri+r2^ p10  *  r2P14 

11 

-(A1+A2)pu  ,  y2P12  ,  pjP13 

12 

A2P11  ’  " (Ai+y2+r2^P12  5  ylP14  5  r2P16 

13 

rlP9  *  A1P11  *  “(A2+y i+ri^pi3  5  y2P14 

14 

rlp10  *  A1P12  *  A2P13  *  “  (y i+y 2+rl+r>2  ^ P14  *  r2P18 

15 

-(a1+A2)p15  ,  y2P16  ,  yjP17 

16 

A2P15  ’  -( Ai+y2+r2 ^P16  5  ylP18 
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TABLE  2  (Cont) 


STATE  EQUATIONS  FOR  A  SIMPLE  SYSTEM 
WITH  3-SPACE  BUFFER 


(See  Table  1  for  definition  of  states). 


Derivative  of 
PIState  No.} 

Terms  in  Right-Hand-Side  Sum 

17 

rlp13  ’  Xlp15  *  "(Vyl+rl)p17  •  y2p18 

18 

rlp14  *  Xlp16  5  X2P17  *  "W+y2+rl+r2^p18  »  r2p20 

19 

rlp17’  "X2P19  ’  y2p20 

20 

r 1P18  *  X2P19  *  ”^y2+r2^p20 

12 
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DRSMC-SAS  (R) 

SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 

18.  The  last  5  system  states,  which  contain  waiting  states  for  operation  1, 
must  be  written  explicitly.  For  the  general  case,  these  are: 


• 

pns-4 

A2pns-5"^l+u2+r2^pns-4 

•Vns-2 

(9) 

Pns-3 

“  rlpns-7+:‘lPNs-5"*:2+ul+rl)pn5-3 

+u2pns-2 

(10) 

pns-2 

'  rlpns-6+'lpns-4+i2pns-3 

-(Ul+y2+r1+r2)pns_2+r2pns 

(11) 

pns-l 

rlpns-3”2pns-  l+!J2pns 

(12) 

pns  = 

rlpns-2n2pms-r(p2+r2)pns  ' 

•  (13) 

state 

transition  diagram  for  the  last  five  states  is  found 

in  Figure  3. 

19.  To  evaluate  the  system  under  stochastic  steady  state,  the  £  vector  is 
set  to  zero. 

Then, 

B£  =  0  .  (14) 

As  noted  above,  the  resultant  set  of  equations  contains  one  superfluous 
equation,  since  the  state  probabilities  must  sum  to  1.  To  remedy  this 
situation,  (14)  is  converted  to  the  linear  matrix-vector  equation 

A  x  =  c_ 

with 

A  (nxn),  x  (nxl),  and  c  (nxl)  , 

where 

..n  =  ns-1  . 

This  is  done  by  assigning 

P1  =  1_Zk=2  pk  ‘  (16) 
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SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 

The  elements  of  the  A  matrix  and  c  vector  in  (15)  are  obtained  via  the 
transformations 

b'ij '  Vbu  •  2i1ins 

and, 

»  2£i£ns 

aij  =  b'i+l,  J+l  ’  ^  -  J i"  • 

(17) 

Note  that  the  first  scalar  equation  in  equation  (14)  is  deleted  in  forming 
(15).  The  A  matrix  is  of  full  rank,  so  £  may  be  obtained  by 

X  =  A~£  . 

(18) 

Then, 

Pi+1  =  xi  ’  l<i<n  • 

(19) 

Finally,  is  obtained  from  (16). 

20.  The  probabilities  of  buffer  occupancy  are  obtained 
probability  vector  £  (or,  alternatively  from  x.)  by 

from  the  state 

Plbuffer  is  empty}  =  p,^  . 

(20a) 

P(j  parts  in  buffer}  = 

Ek=l  pk+4j+2  ,  • 

(20b) 

P(buffer  is  full}  =  (by  definition) 

P(m  Parts  in  buffer}  =  E®=1  Pk+4m+2  . 

(20c) 

The  probability  that  operation  2  must  wait  for  a  buffer 
the  sum  of  the  first  two  state  probabilities: 

condition  is  just 

Ploperation  2  must  wait}  =  p1+p2  . 

(21) 

15 


DRSMC-SAS  (R)  . 

SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 

The  probability  that  operation  1  must  wait  for  a  buffer  condition  is  the  sum 
of  the  last  two  state  probabilities: 

P{operation  1  must  wait}  =  Pn  +  Pn+i  •  ^ 


The  system  productivity  depends  upon  the  probability  that  operation 
down.  This  last  probability  is  given  by 

(ns-2)/2 


Ploperation  2  is  down}  =  E 

k  = 


p2k+l 


2  is 


(23) 


Finally,  the  system  productivity  is 

1  -  P{operation  2  is  down}  -Pfoperation  2  must  wait}  .  (24) 

If  one  wants  the  steady-state  production  rate  of  the  simple  system,  this  value 
is  obtained  from  the  productivity  by  multiplying  times 


21.  Analytic  Results  of  GS.BUF 

Calculations  were  made  using  GS.BUF  to  compare  results  with  stochastic 
simulation  and  to  do  a  sensitivity  analysis  of  certain  parameters.  This 
computational  experience  provided  timing  estimates  on  our  PRIME  550  mini¬ 
computer.  All  calculations  were  made  using  double-precision  arithmetic. 

The  stochastic  simulation  used  for  comparison  was  a  very  simple  implementa¬ 
tion  of  TANDEMT.  In  all  cases  studied  the  time  between  failures  and  the 
time  to  repair  were  exponential  random  variables.  Simulation  runs  were 
40  24-hour  days,  starting  with  an  empty  system.  Runs  were  made  for  instances 
in  which  the  machine  operating  (or  service)  time  is  exponential  and  in  which 
it  is  constant.  The  runs  with  constant  service  time  are  used  to  test  the 
applicability  of  the  Markov  model  in  GS.BUF  to  a  quite  different  model. 
Parenthetically,  I  note  that  other  Markov  models  can  be  created  to  approximate 
the  constant-rate  case.  By  defining  many  substates  to  describe  a  machine, 
operation,  it  is  possible  to  describe  a  random  service  time  whose  coefficient 
of  variation  is  quite  small  (if  not  zero).  In  fact,  I  constructed  a  specific 
case  of  such  a  model  using  3  substates.  The  productivity  calcul ated  with  that 
model  is  in  better  agreement  with  simulated  results  than  is  the  case  for  GS.BUF. 
However,  due  to  the  large  number  of  system  states  produced  by  this  method, 
computational  efficiency  is  poor,  for  a  typical  buffer  size.  Consequently, 
that  approach  was  abandoned.  It  is  recommended  for  single-machine  operations 
having  constant  rates,  that  stochastic  simulation  be  used  to  estimate  system 
productivity.  This  approach  also  has  the  advantage  of  modeling  other-than 
well  maintained  systems,  which  were  considered  here. 
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SUBJECT:  Methods  for  Choosing  Buffer  Size  in  Tandem  Production  Operations 

22.  Results  of  GS.BUF  in  which  the  operations  have  a  common  set  of  parameters 
are  shown  on  the  first  page  of  Table  3.  Buffer  capacity  is  treated  as  a 
parameter  in  the  comparison  between  calculated  and  simulated  results.  This 
comparison  indicates  complete  agreement  within  expected  statistical  variation. 
For  simulations  of  this  length  the  estimated  standard  error  of  the  productivity 
estimate  is  about  0.11  to  0.14  for  these  examples.  The  typically  somewhat 
larger  estimated  probability  that  the  buffer  is  empty  is  a  reflection  that 

the  simulation  was  not  started  in  the  steady  state.  However,  due  to  the 
simulation  length,  this  effect  is  small.  On  the  second  page  of  Table  3,  a 
comparison  is  made  in  which  the  parameters  of  operation  1  are  not  the  same 
as  those  of  operation  2.  Again,  agreement  is  excellent.  Whenever  one  is 
designing  a  buffer,  a  balanced  pair  of  operations  should  be  considered. 

When  balanced,  the  operation  thruput  is  a  constant  equal  to  the  product  of 
the  machine  rate,  the  machine  availability,  and  the  number  of  machines. 

Note  that  these  operations  are  balanced. 

23.  Calculated  and  simulated  productivities  for  an  expanded  set  of  buffer 
capacities  are  shown  in  Table  4.  This  sort  of  analysis  can  be  used  in 
making  economic  tradeoffs  when  choosing  a  buffer  capacity.  The  results 

of  TANDEMT  are  shown  in  Table  5.  It  is  noted  that  the  system  productivity 
is  always  greater  when  the  machine  service  (operating)  times  are  constants 
than  when  they  are  exponential  random  variables.  However,  note  that  the 
difference  (and  relative  difference)  in  productivities  in  these  two  instances 
diminishes  as  the  buffer  capacity  increases.  This  fact  suggests  that  regard¬ 
less  of  the  distribution  of  service  times,  GS.BUF  may  be  a  practical  procedure 
to  use  for  productivity  estimation  (and  tradeoff)  when  the  buffer  capacity 
is  large.  Figure  4  shows  the  probability  density  functions  of  buffer  state 
occupancy  for  cases  in  which  the  service  times  are  constant  and  exponential 
random  variables.  The  U-shaped  densities  are  typical.  Note  that  much  larger 
probabilities  of  being  at  the  buffer  extremes  is  exhibited  by  a  constant-rate 
system. 

24.  The  sensitivity  analysis  using  GS.BUF  considers  the  effect  of  the 
following  three  parameters  on  system  productivity:  buffer  capacity,  MTTR, 
and  machine  rate.  While  not  large,  the  ranges  of  these  parameters  are 
representative  of  many  ammunition  production  operations.  Several  inferences 
can  be  drawn  from  the  results  shown  in  Table  6.  Only  one  is  mentioned  at 
this  point.  Suppose  the  machine  rate  is  given  and  the  buffer  size  chosen 

to  yield  a  particular  productivity  or,  alternatively,  chosen  so  that  the 
marginal  cost  of  additional  buffer  space  just  equals  the  value  of  additional 
production  from  a  system  with  the  incrementally  larger  buffer.  If,  later, 
a  greater  machine  rate  is  available  via,  possibly,  machine  substitution,  one 
can  expect  a  productivity  decrease  if  the  buffer  is  not  enlarged.  Remember 
that  productivity  is  a  measure  of  production  efficiency.  Thus,  doubling  the 
machine  rate  —  with  buffer  fixed  —  will  increase  the  production  rate,  but 
will  not  double  it. 
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TABLE  3 


COMPARISON  OF  THE  PROBABILITIES  OF  BUFFER 
OCCUPANCY:  CALCULATED  VERSUS  SIMULATED 


Parameters: 

Two  single-machine  operations  in  tandem. 
Operation  (service)  times  are  exponential. 
Repair  times  are  exponential. 

MTBF  =  100  minutes 
MTTR  =  25  minutes 

Average  operating  rate  =  1  part/min. 


Buffer 

Buffer 

State 

Probabil ities 

Capacity 

Status 

Calcul a ted 

Simulated 

3 

0 

0.379 

0.385 

1 

0.121 

0.125 

2 

0.121 

0.126 

3 

0.379 

0.364 

System  Productivity* 

0.595 

0.599 

5 

0 

0.320 

0.312 

1 

0.091 

0.094 

2 

0.089 

0.095 

3 

0.089 

0.094 

4 

0.091 

0.094 

5 

0.320 

0.311 

System  Productivity* 

0.620 

0.630 

10 

0 

0.250 

0.247 

1 

0.059 

0.061 

2 

0.056 

0.059 

3 

0.055 

0.056 

4 

0.054 

0.054 

5 

0.053 

0.056 

6 

0.054 

0.054 

7 

0.055 

0.056 

8 

0.056 

0.057 

9 

0.059 

0.060 

10 

0.250 

0.240 

System  Productivity* 

0.650 

0.640 
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TABLE  3  (Cont) 


COMPARISON  OF  THE  PROBABILITIES  OF  BUFFER 
OCCUPANCY:  CALCULATED  VERSUS  SIMULATED 


Parameters: 

Two  single-machine  operations  in  tandem. 
Operation  (service)  times  are  exponential. 
Repair  times  are  exponential. 

MTBF  =  200  min  (1st  opn),  =  100  min  (2nd  opn). 
MTTR  =  25  min  (1st  opn),  =  12.5  min  (2nd  opn). 
Average  operating  rate  =  1/min. 


Buffer 

Capacity 

Buffer 

States 

State 

Cal culated 

Probabil ities 

Simulated 

3 

0 

0.357 

0.368 

1 

0.140 

0.143 

2 

0.141 

0.144 

3 

0.362 

0.345 

System  Productivity* 

0.695 

0.696 

5 

0 

0.287 

0.291 

1 

0.103 

0.105 

2 

0.104 

0.106 

3 

0.105 

0.106 

4 

0.108 

0.109 

5 

0.294 

0.283 

System  Productivity* 

0.727 

0.731 

*  System  productivity  is  the  ratio  of  the  average  production  achieved  to 
the  maximum  steady-state  production  from  a  system  of  perfect  machines 
operating  at  the  same  rates. 
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TABLE  4 


COMPARISON  OF  CALCULATED  WITH 
SIMULATED  PRODUCTIVITY  ESTIMATES  FOR 
A  SIMPLE  PRODUCTION  SYSTEM* 


Buffer 

Average  Productivity 

Capacity 

Calculated 

Simul ated** 

3 

0.595 

0.599 

5 

0.620 

0.630 

10 

0.650 

0.640 

20 

0.680 

0.684 

40 

0.710 

0.713 

*  Parameters: 

A  single  machine  at  each  of  two  operations. 

Exponential  operation  times. 

Exponential  repair  times. 

Average  operating  rate  1  part/minute. 

Common  MTBF  =  100  minutes. 

Common  MtTR  =  25  minutes. 

**  The  standard  error  is  about  0.011,  based  on  a  40  day  simulation. 
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TABLE  5 


COMPARISON  OF  PRODUCTIVITY  ESTIMATES  FOR 
SIMPLE  PRODUCTION  SYSTEMS*  HAVING 
CONSTANT  VERSUS  EXPONENTIAL  OPERATING  TIMES 


Tabulated  productivity  is  the  ratio  of  expected  production  to  the  maximum 
production  from  a  perfect  system  operating  at  the  same  machine  rate. 


Buffer 

Capacity 

Simulated  Operating 

Constant 

Times 

Exponential 

3 

0.674 

0.599 

5 

0.679 

0.630 

10 

0.690 

0.640 

20 

0.709 

0.684 

40 

0.724 

0.713 

*  Parameters: 

Single  machine  operations. 
Average  operating  rate  1/min. 
Repair  times  are  exponential. 
MTBF  =  100  minutes 
MTTR  =  25  minutes 
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TABLE  6 


CALCULATED  PRODUCTIVITY  OF  A  SIMPLE  PRODUCTION 
SYSTEM*  AS  A  FUNCTION  OF  SEVERAL  VARIABLES 


The  tabulated  productivity  is  the  ratio  of  average  production  achieved  to 
the  maximum  steady-state  production  from  a  perfect  system  operating  at  the 
same  machine  rate. 


Buffer 

Capacity 

MTTR 

(minutes) 

Machine  Rate  (parts/min) 

1 

2 

10 

12.5 

0.770 

0.761 

25.0 

0.650 

0.642 

20 

12.5 

0.805 

0.792 

25.0 

0.680 

0.666 

40 

12.5 

0.835 

0.819 

25.0 

0.710 

0.690 

*  System  parameters: 

Two  single-machine  operations  in  tandem. 
Capacity  of  intermediate  buffer  is  a  parameter. 
Common  machine  rate  is  a  parameter. 

Common  MTTR  is  a  parameter. 

Exponential  operation  (service)  times. 
Exponential  repair  times. 

Common  MTBF  =  100  minutes. 
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25.  The  system  model  of  GS.BUF  assumes  single-machine  operations.  However, 
multi-machine  operations  can  be  approximated  by  this  model  by  scaling  the 
machine  rate  to  Nr,  where  N  is  the  number  of  machines  in  the  operation, 

each  having  rate  r.  This  approximation  exploits  a  theorem  of  random  processes 
on  the  pooling  of  many  component  processes.  Cox*  states  that  when  many 
independent  process  events  are  pooled,  the  pooled  process  is  approximately 
Poisson,  i.e.,  the  time  between  events  is  approximately  exponential.  This 
is  true  irrespective  of  the  distributions  of  the  component  processes.  In 
application  to  multi-machine  operations,  each  machine's  output  gets  pooled 
for  that  operation.  Thus,  the  times  between  unit  production  events  are 
approximately  exponential  random  variables  when  the  number  of  machines  is 
large.  Surprisingly,  the  number  working  in  parallel  at  an  operation  does 
not  need  to  be  more  than  about  4  to  yield  a  good  approximation  of  the 
probability  density  function  (p.d.f.)  for  buffer  occupancy.  This  point 
is  illustrated  by  the  results  in  Table  7.  Two  multi-machine  cases  were 
simulated:  4  machines  per  operation  working  at  machine  rate  1/4  parts 
per  minute  and  5  machines  working  in  parallel  at  machine  rate  1/5  parts 
per  minute.  For  comparison  is  shown  the  p.d.f.  of  buffer  occupancy 
calculated  with  GS.BUF  having  a  machine  rate  of  1  part  per  minute.  One 
observes  that  the  p.d.f.'s  in  these  instances  are  nearly  the  same. 

Agreement  between  calculated  and  simulated  productivities  is  not  as  good, 
however.  The  simulated  productivities  for  the  two  examples  are  both  0.87 
whereas  the  calculated  productivity  is  0.81. 

26.  Conclusions  Regarding  GS.BUF 

Specific  conclusions  for  the  method  of  GS.BUF  are  summarized  here. 

(a)  When  a  simple  production  system  satisfies  the  assumption  of  exponentially 
distributed  service  times,  the  results  of  GS.BUF  are  exact.  This  point  has 
been  verified  with  simulation,  (b)  With  stochastic  steady  state,  the 
probability  density  function  (p.d.f.)  of  buffer  occupancy  is  U-shaped.  For 
operations  having  a  common  set  of  operating  characteristics,  the  above  p.d.f. 
is  symmetric  with  respect  to  0.5 -buffer  capacity  spaces.  A  prominent 
positive  jump  in  the  p.d.f.  is  observed  at  the  first  and  last  states.  For 
a  given  buffer  capacity  the  jump  is  more  pronounced  when  the  operating  times 
are  constant  than  when  they  are  both  exponential.  This  implies  that 
relatively  more  time  is  spent  at  extreme  states  when  the  operating  times 
are  constant  than  when  they  are  random.  Nevertheless,  the  productivity  of 
a  system  with  constant  operating  rate  is  greater  than  that  of  a  system  with 
a  random  rate  of  the  same  mean  value,  other  things  being  the  same. 


*  Page  77  ff.  Cox,  D.R.  Renewal  Theory,  London,  distributed  by  Barnes 
and  Nob! e,  c .  1962. 
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TABLE  7 


COMPARISON  OF  AN  ANALYTIC  APPROXIMATION*  WITH  SIMULATED 
PROBABILITY  DISTRIBUTIONS  OF  BUFFER  OCCUPANCY  FOR 
SYSTEMS  HAVING  MULTIPLE  MACHINES  PER  OPERATION 


*  Parameters  of  the  Analytic  Model : 

Two  single-machine  operations  in  tandem. 
Operation  times  are  exponential. 

Repair  times  are  exponential. 

MTBF  for  each  machine  =  100  minutes. 

MTTR  for  each  machine  =  12.5  minutes. 
Average  operating  rate  =  1  part/mi n. 
Buffer  capacity  =  20  spaces. 


Buffer 

State 

Cal  c. 
pdf 

-f* 

Simulated 

4  mach  at  rate  1/4 

pdf  with 

5  mach  at  rate  1/5 

0 

0.141 

0.131 

0.135 

1 

0.044 

0.037 

0.042 

2 

0.042 

0.035 

0.043 

3 

0.040 

0.036 

0.044 

4 

0.038 

0.037 

0.043 

5 

0.037 

0.037 

0.041 

6 

0.036 

0.041 

0.036 

7 

0.036 

0.042 

0.036 

8 

0.035 

0.041 

0.038 

9 

0.035 

0.041 

0.036 

10 

0.035 

0.039 

0.033 

11 

0.035 

0.035 

0.034 

12 

0.035 

0.039 

0.034 

13 

0.036 

0.040 

0.035 

14 

0.036 

0.040 

0.037 

15 

0.037 

0.040 

0.039 

16 

0.038 

0.037 

0.037 

17 

0.040 

0.035 

0.036 

18 

0.042 

0.037 

0.040 

19 

0.044 

0.037 

0.043 

20 

0.141 

0.143 

0.138 

+  Simulated  machine  operating  rates  are  constants.  The  MTBF  and  MTTR 
parameters  of  the  simulation  are  as  given  above. 
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(c)  A  sensitivity  analysis  of  system  parameters  was  performed  using  GS.BUF. 
All  the  following  are  shown  to  be  important:  Machine  operating  rate  (or 
service  time),  repair  time,  and  buffer  capacity,  (d)  If  the  buffer  capacity 
remains  constant  and  the  machine  rate  is  doubled,  a  productivity  decrease  is 
experienced  --  even  tho  the  production  rate  of  the  new  configuration  is 
greater.  This  loss  of  productivity  is  more  pronounced  for  a  system  with  a 
large  buffer  than  with  a  system  whose  buffer  is  inadequate,  (e)  Doubling 
the  production  rate  and  concurrently  doubling  buffer  capacity  increases 
the  system  productivity.  Thus,  the  buffer  size  need  not  be  doubled  to 
preserve  productivity  if  the  operating  rate  is  doubled,  (f)  At  any  reason¬ 
able  buffer  size,  halving  the  MTTR  produces  a  greater  improvement  in 
productivity  than  that  achieved  by  doubling  buffer  capacity,  (g)  For 
systems  having  constant  machine  operating  rates,  GS.BUF  can  be  used  to 
approximate  the  system,  only  if  many  machines  are  working  in  parallel  at 
each  operation.  The  probability  distribution  of  buffer  state  occupancy  is 
approximately  correct  in  this  case.  If  only  one  fixed-rate  machine  exists 
per  operation,  it  is  recommended  that  stochastic  simulation  be  used  to 
estimate  the  system  productivity. 

27.  Derivation  of  Equations  for  BUF.CAP 

To  motivate  subsequent  discussion,  consider  the  state  dynamics  of  a  one- 
machine  system.  This  system  is  regarded  as  operating  independently  of  other 
buffers  and  production  operations.  In  this  case  there  are  only  two  states  -- 
down  (0)  and  up  (1).  Markov  transitions  from  down  to  up  occur  at  the  birth 
rate  X,  and  transitions  from  up  to  down  occur  at  the  death  rate  y.  Thus, 
the  equations  for  the  probabilities  of  state  occupancy  are 

P0(t)  =  -XPq+mPj  (25a) 

Pjtt)  =  XpQ-yp1  .  (25b) 


These  equations  can  be  solved  directly  or  by  using  Laplace 
i nverting. 

P0(t)  *  p0(0)e_(x+u)t  +  (1-e'(Xt',)t) 

Pl(t)  =  (l-e-(x+*,)t)  +  p1(0)e-<x+*l)t  . 

Note  that  the  steady-state  availability, 

P1(~) >  or 


A  = 


X 

X+y 


transforms  and 

(26a) 

(26b) 

(27) 
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Using  the  definitions  of  A  and  u,  equations  (1,  2), 


MTBF 

MTBF+MTTR 


(28) 


The  (mathemati sally)  expected  production  from  this  operation  over  the  time 
interval  (0,  t)  is  given  by 


Pj(x)dx 


(29) 


where  r  is  the  machine  operating  rate. 

From  (26b)  and  (29), 

-  _2_  t  .  X  , 

r  A+y  t  (a+u)2 


P.i(Q) 

A+y 


)  • 


(30) 


28.  Consider  two  single-machine  operations  in  series  with  a  common  set  of 
operational  parameters.  When  first  observed,  let  the  first  operation  be  in 
state  1  and  the  second  operation  be  in  state  2.  Call  this  initial  condition 
IC1.  Stated  mathematically,  let 

Pj(0)  =  1  for  operation  1 
and 

Pj(0)  =  0  for  operation  2. 

Then,  the  difference  in  expected  production  of  these  operations  as  t  becomes 
large  can  be  obtained  from  (30)  as  E[production  difference,  given  IC1]  =  r/(A+u). 
(Notationally ,  E  is  the  expected  value  operator.)  If  IC1  obtains,  one  expects  to 
add  the  above  production  quantity  to  the  contents  of  a  buffer  between  the 
operations,  given  the  assumed  freedom  from  buffer  constraints.  Clearly,  if  the 
initial  operation  states  had  been  the  l's  complements  of  the  above  (IC2),  a 
quantity  of  production  would  have  been  removed  from  the  buffer  equal  to  (31). 

This  again  assumes  that  at  least  that  much  was  present  initially.  The  IC's 
considered  are  the  extreme  conditions  for  this  example.  Therefore,  at  reasonably 
low  risk  of  buffer  inadequacy  one  might  suggest  using  the  range  of  E[production 
difference]  over  extreme  IC's  for  sizing  the  buffer.  In  this  case 

Etproduction  difference,  ICl]-E[production  difference,  IC2] 

=  2r/(A+y)  .  (32) 


(31) 
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However,  this  approach  would  ignore  some  obvious  facts. _  First,  we  are 
dealing  with  conditional  expected  values.  Thus,  variation  from  these 
averages  must  occur.  Secondly,  we  ignore  the  fact  that  the  occupancy  of 
the  buffer  at  the  commencement  of  delivery  of  additional  product  under 
I  Cl  may  exceed  r/(A+y).  To  reduce  the  risk  associated  with  these 
contingencies,  one  can  identify  the  variation  in  expected  buffer  state 
change.  Specifically,  one  can  calculate  the  standard  deviation  (SD)  of 
this  random  variable.  Then,  some  number  of  SD's  can  be  added  to  the 
range  in  (32)  to  produce  a  buffer  capacity  requirement.  I  have  somewhat 
arbitrarily  chosen  4  SD's  as  a  reasonably  cautious  number.  The 
probability  distribution  of  expected  production  difference  is  calculated 
by  enumerating  all  possible  IC's.  In  this  example  there  are  four  — 

(0,  0),  (0,  1),  (1,  0),  and  (1,  1).  The  probability  that  the  system 
of  two  operations  exists  in  one  of  these  states  is  simply  the  product 
of  the  probability  that  operation  1  is  in  its  given  state  and  the 
probability  that  operation  2  is  in  its  state.  The  probability  that 
a  one-machine  operation  is  up  (1)  is  A,  and  down  (0)  is  (1-A).  In 
general,  in  a  N-machine  operation 

P{k  machines  are  operating} 

=  (]J)  Ak  (l-A)N"k  .  (33) 

29.  Clearly,  the  expected  production  difference  between  operations  started 
in  the  same  state  is  zero.  Hence,  in  this  example  the  (0,  0)  and  (1,  1)  IC's 
yield  E[production  difference]  =  0.  Thus,  the  probability  density  of  the 
E [production  difference]  for  this  example  is  the  following:  -r/(A+y)  with 
probability  A(l-A),  0  with  probability  A2+(1-A)s  and  r/(A+y)  with  probability 
A(l-A).  Since  the  mean  value  of  this  variable  is  zero,  the  variance  is  given 
as 

[2r2/(A+y)2]A(l-A)  (34) 

or,  from  (27),  the  variance  of  E[production  difference]  is 

2Ay  ^  (35) 

( A+y) 

and  the  SD  is 

(2Xp)1/2r  .  (36) 

U+p)2 

As  a  specific  numerical  example,  let  r=l  part  per  minute,  MTBF  =  100  minutes. 
MTTR  =  25  minutes.  Then,  from  (32),  the  expected  range  in  production 
differences  is  40  parts  and,  from  (36),  the  standard  deviation  11.3  parts. 
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This  yields  a  required  capacity  of  85  spaces.  To  examine  the  marginal 
advantage  in  productivity  for  this  buffer  size,  I  ran  GS.BUF  at  a  buffer 
capacity  of  85  and,  again  at  80.  The  productivities  in  these  instances 
are  0.7426  and  0.7402,  respectively.  This  implies  a  productivity  change 
of  about  0  04%  per  percent  change  in  buffer  capacity. 

30.  A  Two-Machine  Operation 

As  in  the  above  example,  it  is  not  difficult  to  obtain  an  analytic 
solution  to  the  time-dependent  Markov  model  of  a  2-machine  production 
system.  With  two  machines  the  system  states  are  0,  1,  and  2.  However, 
since 

pQ(t)  =  1  -  Pj(t)  -P2(t)  ,  (37) 

only  two  equations  are  required  to  describe  the  system.  Dropping  the 
notation  for  explicit  time  dependence, 

Pj_  =  2A(l-p1-p2)  -  (x+pjpj  +  2yp2  (38a) 

p2  =  Xpj  -  2yp2  .  (38b) 

After  some  manipulation  of  these  equations  and  using  Laplace  transforms. 


one  obtains 

p  (t)  =  lAl*  +  Ae'(A+y)t  +  Be"2(A+y)t  , 

1  .  (A+y)2 

(39a) 

with 

A  _  ao_ai(A+^  _  4 Ay 

A+y  (A+y)2 

(39b) 

Q  _  aQ-2a1(A+y)  +  2Xy 

A+y  (A+y)2 

(39c) 

aQ  =  Pj(0)  +  3(A+y)p1(0) 

or 

aQ  =  2ApQ(0)  +  2(A+y)p1(0)  +  2yp2(0) 

(39d) 

a1  =  P^O)  . 

(39e) 
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2 

Po(t)  =  Po(0)e"2yt+  2 

L  (A+y) 


Buffer  Size  in  Tandem 


A2  e-2yt 
( A+y ) 


+ 


Production  Operations 


A  (e-2yt  _  g-(A+y)t 
A-y  k 

Steady-state  results  are 
approach  °°. 

2 


P0H 

P^00) 

P2(°°) 


y 

(A+y)2 

2Ay 

(A+y)2 

>2 

(A+y)2 


)  +  |  (e-2yt  _  e-2(A+y)t}  .  (40) 

obtained  from  (37),  (39),  and  (40)  by  allowing  t  to 

(41a) 

(41b) 

,  (41c) 


The  expected  number  of  machines  operating  in  steady-state  in  this  case  is 

E[number  operating]  =  p^00)  +  2p2(°°)  »  (42a) 

=  2A/ (A+y)  .  (42b) 


Since  the  number  of  machines,  N,  is  2  in  this  example,  as  anticipated, 

E[number  operating]  =  NA  .  (43) 

The  expected  number  of  machines  operating  equals  the  max  number  times  the 
intrinsic  availability  only  i f  the  system  is  well  maintained.  In  writing 
the  system  equations  (38a  and  b),  it  is  assumed  that  if  N  machines  are 
down  for  repair,  N  repairmen  are  working.  Thus,  the  transition  rate  from 
the  0  state  to  the  1  state  is  given  as  NA.  If  there  were  only  M  repairmen 
available  for  machine  repairs,  where  M<N,  the  largest  transition  rate  from 
a  lower  to  a  higher  state  would  be  MA.  In  the  latter  case,  the  system 
state  probabilities  would  depend  upon  both  the  number  of  machines  per 
operation  and  the  number  of  repairmen.  To  avoid  this  complication,  it  is 
assumed  that  the  system  is  well  maintained.  This  assumption  is  not  as 
restrictive  as  it  may  seem.  Approximately  correct  results  are  obtained 
for  far  weaker  conditions,  as  will  be  shown  later. 
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31.  An  N-Machine  Operation 


m_j. _  r^e  state  transition  diagram  for  the  general  case  is  shown  in  Figure  5. 

rate  N  f  ?  !?°m  the  j  bh  t0  the  j+1  st  states  occurs  at  the 

tho  ir  i  +  +neS  down^  tim®s  The  downward  transition  from  the  k  th  to 

ratP  „  VlatefCCiJr/  &luthe  rat-e  k  (more  t0  fai>1)  times  the  unit  death 

t  '  AS  waf  done  for  the  previous  Markov  models,  the  Kolmogorov 

thlS  can  be  wri’tten  by  inspection  from  the  state 

not  rJnl^p^h^3"1*  tV"™  th\s  Process  is  straightforward,  it  is 

K  nnfn  J  ?  Jhe*ge,\e,:al  result>  with  the  deletion  of  the  zero 
(or  null)  state  is  the  familiar  form 


£(t)  ■  A£(t)  +  £  , 
with  c'  =  [NX,  0,  0 . 0]  , 


(44) 


where  £  and  c  are  (Nxl)  and  where  A  is  (NxN). 

rule  withStime1steprh^edUre  employs  numerical  Integration  using  the  rectangle 


£(t+h)  =  £(t)  +  h£(t) 


(45) 


trhP^iLhtiT^J44)’  .F°r  a  s?a11  time  steP  (h)’  this  Procedure  was  found 
taster  ^a"Q  easler)  t0  implement  than  to  use  a  double  step 
'  with  Euler  s  rule  with  a  predictor  and  a  corrector.  A  time  step  h  of  0  1 
minute  is  used  in  BUF.CAP.  Notationally,  let  the  conditional  expected  value’ 
of  the  output  from  an  operation  at  time  t  (from  the  IC)  be  denoted  7(t). 
by  definition, 


%  ft  N 

x(t)  =  rz  kp.  (t)dt  , 
J ok=i  K 


(45) 


Xfe  Vs  the  unit  machine  rate  for  this  operation.  The  integrand  in  (45)  is 
^•f^-age-rate  °!  Pr°doction  -  dx/dt  -  from  this  machine  operation.  This 
derivative  is  saved  at  time  intervals  of  a  for  optional  printing.  Call  the 

intpnr^t!nap?r0x1!1ati10n  with  inte9er  i.  In  performing  the  numerical 

integration  to  wlculate  x(t),  it  is  unnecessary  to  use  a  step  as  small  as 

n.  T°  yield  about  the  same  relative  precision  as  obtained  in  calculating 
£  t),  one  can  use  a  step  a  -  called  DELTAT  in  BUF.CAP  -  of  0.5  minute 
with  Euler  s  rule: 


X ( t  +A )  =  0.5A[X(t+A)  +  £(t)l 


(46) 
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Denote  by  x-j-j  the  conditional  expected  output  of  operation  i  with  initial 
condition  J  (IC)  j.  Then,  the  difference  in  expected  production  values 
is  calculated  at  time  t  : 

ITla  A 


Axj  "/lj^max^  "  x2j^max^ 


(47) 


This  difference  is  essentially  constant  beyond  tmax  =  2max  (MTTRj,  MTTRg). 
These  expected  production  differences  are  treated  in  BUF.CAP  in  the 
manner  described  for  the  single-machine  case  (pgf.  28  ff).  The 
probability  distribution  of  TTxj  is  displayed  in  order  to  facilitate  the. 
choice  of  buffer  capacity  on  tne  basis  of  risk  that  a  particular  value  is 
inadequate. 


33.  Sample  Output  from  BUF.CAP 

A  sample  run  using  the  program  BUF.CAP  is  shown  in  Table  8.  The  program 
input  values  are  repeated  at  the  top  of  the  page.  The  output  represents  an 
abbreviated  form  of  the  available  outputs.  If  the  user  chooses  to  display 
the  trajectories  of  the  conditional  state  probability  vectors  and  cumulative 
production  he  may. 

34.  Stochastic  simulation  is  used  in  the  process  of  examining  the  advantage 
of  increasing  the  buffer  capacity  beyond  that  required  by  BUF.CAP.  Examples 
of  the  c.d.f. 's  of  buffer  occupancy  for  several  two-operation  systems  are 
shown  in  Figures  6  and  7.  The  max  operation  rate  of  both  operations  in  each 
system  is  the  same.  In  all  cases  shown  here  the  machine  service  (or  operation) 
time  is  constant.  The  probability  distributions  of  buffer  occupancy  for  two 
balanced  systems  are  compared  in  Figure  6.  One  can  observe  that  increasing 

the  buffer  capacity  from  40  to  100  spaces,  for  the  parameters  shown,  has  the 
effect  of  significantly  reducing  the  risk  of  encountering  a  full  buffer.  But, 
the  probability  of  an  empty  buffer  is  nearly  the  same  in  both  instances. 

Several  comparisons  are  made  in  Figure  7.  In  these  comparisons  the  buffer 
capacity  is  fixed  at  40  spaces,  and  the  operational  rate  is  a  constant  1  part/ 
minute.  An  effect  on  the  c.d.f.  of  buffer  occupancy  occurs  when  the  number  of 
machines  per  operation  increases.  The  effect  of  this  increase  is  to  reduce 
the  probabilities  associated  with  the  extreme  states  of  the  buffer.  This 
phenomenon  is  also  shown  in  the  output  of  BUF.CAP.  Another  observation  of 
interest  can  be  drawn  from  Figure  7.  Thruout  this  study  all  analyses  have 
been  conducted  assuming  a  well  maintained  system.  As  indicated,  this  implies 
at  least  as  many  repairmen  as  machines.  For  comparison,  one  simulation  run 
was  made  with  a  system  of  6  machines  —  4,  in  operation  1  and  2,  in  operation  2-- 
and  with  only  one  repairman.  The  c.d.f.  of  buffer  occupancy  for  this  case  is 
shown  in  Figure  7.  A  rather  small  difference  exists  between  the  c.d.f.  for 
this  case  and  the  c.d.f.  for  a  comparable,  well -maintained  system.  The  reason 
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TABLE  8 


SAMPLE  OUTPUT  FROM  PROGRAM  HUP. CAP  WITH  TERMINAL  DIALOG  DELETED 


machine  input  data  for  buffer  capacity  calculation 

OPERATION  1  OPERATION  2 


FIE  machines  _ 

4 

4 

MACH  RATE 

0  .2  5 

0 

.25 

PARTS/MIN 

Ml  B  F 

i  o  o .  o  e 

100 

.0  0 

MINUTES 

M  T  T  K  ‘ 

li  *50 

12 

.50 

MINUTES 

availability 

88.89 

88 

.89 

PERCENT 

feptht 

"TThE”FTnT5TES7_'3!)7 

u - *ssormro 

~STFT_EWmENEE  q7^2  ” 

EXPECTED  BUFFER  REQUIREMENTS  ORDERED  OVER  ALL  SYSTEM  STATES 

BASED  ON  A 

REPAIR  TIME 

LAG  OF 

50.0  MINUTES 

CTITBeR 

FufFFr" 

— utfeer — 

TFoE - 

CUM  L 

- OP  N  T 

~“"5Fn  1 

INDEX 

SPACES 

SPACES 

DENS 

PkOb 

STATE 

STATE 

- T - 

”-T77F" 

- 77Ti - 

■(T777IFT“ 

HrannrT- 

TJ 

4 

2 

-8.2 

2.7 

0.0030 

0  .0  031 

1 

4 

3 

-8.2 

2.7 

0.0000 

0.0032 

0 

3 

4 

-5.5 

5.5 

0.0000 

0.0032 

0 

2 

5 

-5.5 

5.5 

0.0365 

0.0397 

2 

4 

6 

-5.5 

5.5 

0.0015 

0.0413 

1 

3 

7 

-2.7 

8.2 

0.1949 

0  .2361 

3 

4 

8 

-2.7 

8.2 

0.0000 

0.2361 

0 

1 

9 

-2.7 

8.2 

0.0003 

0  .2364 

1 

2 

10 

-2.7 

8.2 

0.0183 

0  .2547 

2 

3 

11 

0.0 

10.9 

0.0000 

0  .2547 

1 

1 

12 

0.0 

10.9 

0.0000 

0  .2547 

0 

0 

13 

0.0 

10.9 

0.0974 

0  .3521 

3 

3 

14 

0.0 

10.9 

0.0034 

0  .3556 

2 

2 

15 

0.0 

10.9 

0.3897 

0.7453 

4 

4 

16 

2.7 

13.7 

0.0183 

0  .7636 

3 

2 

17 

2.7 

13.7 

0.0003 

0  .7639 

2 

1 

18 

2.7 

13.7 

0.0000 

0.7639 

1 

0 

19 

2.7 

13 .7 

0.1949 

0.9587 

4 

3 

20 

5.5 

16.4 

0.0015 

0.9603 

3 

1 

21 

5.5 

16.4 

0 . 0365 

0  .9968 

4 

2 

22 

5.5 

16.4 

0.0000 

0.9968 

2 

0 

23 

8.2 

19.1 

0.0000 

0.9969 

3 

0 

24 

8.2 

19 .1 

0.0030 

0  .9999 

4 

1 

25 

10.9 

21.9 

0.0001 

1.0000 

4 

0 

AVERAGE  BUFFER  CHANGE  -0.00 

STD  DEV  BUFFER  CHANGE  2.43 

REQUIRED  BUFFER  CAPACITY  =  32 
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for  the  slight  effect  of  additional  repairmen  here  is  that  for  this  system 
the  probability  that  other  than  one  repairman  is  needed  is  only  about  12%. 

Thus,  the  assumption  of  the  analysis  that  the  system  is  well  maintained, 
may  not  be  as  restrictive  as  it  first  seems. 

35.  Conclusions  Regarding  BUF.CAP 

(a)  The  risk  of  exceeding  the  buffer  capacity  requirement  calculated 
by  BUF.CAP  is  generally  quite  small  —  typically  less  than  10%.  Further, 
the  marginal  productivity  change,  evaluated  at  the  required  capacity,  is 
nearly  a  constant  0.04%  per  %  change  in  buffer  capacity,  (b)  For  a  large  — 
say,  >80  spaces  --  buffer,  BUF.CAP  executes  faster  than  GS.BUF.  This  may  be 
a  consideration  for  execution  on  small  computers.  Actually,  the  execution 
time  of  BUF.CAP  depends  on  the  number  of  machines  at  each  operation,  not  on 
the  required  buffer  capacity  as  such,  (c)  Unless  the  number  of  repairmen 
is  at  least  equal  to  the  total  number  of  machines,  there  is  a  finite 
probability  that  machines  must  queue  for  repairs.  If  this  happens,  machine 
availability  is  not  equal  to  intrinsic  availability  (A)  and  the  expected 
number  of  machines  is  not  equal  to  N  A.  However,  this  situation  is  not  as 
restrictive  as  it  may  seem.  Both  simulation  and  BUF.CAP  show  that  the 
probability  of  exceeding  a  given  buffer  size  decreases  as  the  number  of 
machines  in  each  operation  increases,  at  constant  thruput.  Under  certain 
conditions  the  probability  distribution  of  buffer  occupancy  is  nearly 
unchanged  by  an  increase  in  the  number  of  repairmen.  This  occurs  at  a  value 
of  number  of  repairmen  M  such  that  prob  (busy  repairmen  <=M)>0.9,  approximately. 
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ANNEXES 


COMPUTER  SOURCE  PROGRAMS 


The  program  listings  in  these  annexes  are  written  in  SIMSCRIPT  2.5.  They 
do  not  employ  any  features  unique  to  the  PRIME  550  minicomputer  on  which  they 
were  run.  Cross-reference  lists  are  included  with  the  program  statements  to 
facilitate  the  identification  of  variable  type  and  locations  within  the 
program.  Since  SIMSCRIPT  is  a  language  very  English-like,  programs  can  be 
followed  easily  without  a  flow  chart.  Therefore,  no  such  diagrams  are 
included.  However,  the  major  program  blocks  are  announced  via  comments, 
which  are  distinguished  from  executable  code  by  stating  with  double  quote 
marks. 

Potential  users  of  these  programs  who  do  not  have  SIMSCRIPT  compilers 
but  do  have  FORTRAN  compilers  are  assured  that  conversion  to  FORTRAN  is 
straightforward.  The  code  in  FORTRAN  is  not  much  longer  than  that  in 
SIMSCRIPT.  If  a  FORTRAN  version  is  implemented  on  a  32  bit  (or  less) 
machine,  it  is  recommended  that  double-precision  arithmetic  be  used.  This 
is  necessary  to  avoid  truncation  error  when  inverting  large  matrices. 

Both  main  programs  were  designed  for  running  interactively.  Inputs  are 
assigned  following  prompting  messages  sent  to  the  terminal. 
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ANNEX  1 


PROGRAM  GS.BUF 


A  sketch  of  the  method  used  to  calculate  the  stochastic  steady  state 
of  a  simple  production  system  is  provided  in  the  body  of  this  memorandum 
under  "Methodology  Overview."  Detailed  system  equations  are  derived  in 
the  section  "Derivation  of  Equations  for  GS.BUF." 

Input  data  is  provided  from  the  terminal  in  response  to  prompting 
messages  such  as  "Input  the  operating  rate  for  the  1st  operation  in  parts 
per  minute."  Output  is  sent  directly  to  the  terminal  for  display.  Since 
this  output  is  often  lengthy,  it  is  recommended  that  a  COMO  file  be 
established  to  display  or  print  it. 
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ANNEX  2 


PROGRAM  BUF.CAP 


BUF.CAP  calculates  a  recommended  capacity  for  a  buffer  separating  two 
tandem  production  operations.  The  methods  employed  in  BUF.CAP  are  outlined 
in  the  body  of  this  memorandum  under  "A  Second  Approach."  The  system 
equations  for  each  production  operation  are  displayed  in  various  places 
depending  on  the  number  of  machines  (N)  in  a  given  operation.  For  N=l, 
see  equation  (26);  for  N=2,  see  (39,  40);  and  for  N>2,  a  general  form  is 
provided  in  equation  (44).  The  elements  of  the  matrix  A  in  (44)  are  given 
in  the  source  program  in  LET  statements. 

Input  data  is  provided  from  the  terminal  in  response  to  prompting 
messages  such  as  "Input  the  number  of  machines  in  1st  operation."  To 
assure  a  system  balance,  the  machine  rates  assumed  in  BUF.CAP,  i.e. 
calculated  endogenously,  may  not  equal  the  actual  or  desired  rates.  TO 
account  for  this  absence  of  input  rates,  the  buffer  requirement  from 
BUF.CAP  must  be  scaled  in  proportion  to  the  ratio  of  actual  thruput  to 
assumed  thruput.  The  program  output  is  sent  to  the  terminal  for  display. 
This  includes  an  echo  of  all  program  inputs.  No  output  files  are  created. 
If  the  output  is  to  be  saved,  a  COMO  file  must  be  established. 
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